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Context. The gas and dust content of protoplanetary disks evolves over the course of a few Myr. As the stellar X-ray and UV light 
penetration of the disk depends sensitively on the dust properties, trace molecular species like HCO^, HCN, and CN are expected to 
show marked differences from photoprocessing effects as the dust content in the disk evolves. 

OAims. We investigate specifically the evolution of the UV irradiation of the molecular gas in protoplanetary disks around a sample of 
' classical T Tauri stars in Taurus that exhibit a wide range in grain growth and dust settling properties. 

Methods. We obtained HCO+ (J = 3-2), HCN (J = 3-2), and CN (J = 2-1) observations of 13 sources with the James Clerk Maxwell 
Telescope. Our sample has 1.3 mm fluxes in excess of 75 mjy, indicating the presence of significant dust reservoirs; a range of dust 
settling as traced through their spectral slopes between 6, 13, and 25 fim; varying degrees of grain growth as extrapolated from the 
strength of their 10-yum silicate emission feature; and complements data from the literature, essentially completing the molecular 
' line coverage for the 21 brightest millimeter targets in the Taurus star-forming region. We compare the emission line strengths with 

, the sources' continuum flux and infrared features, and use detailed modeling based on two different model prescriptions to compare 

\^ ■ typical disk abundances for HCO^, HCN, and CN with the gas-line observations for our sample. 

Results. We detected HCO+ (3-2) toward 6 disks, HCN (3-2) from disks, and CN (2-1) toward 4 disks, with typical 3cr upper limits 
of 150-300 mK (T^f,) in 0.2 km s"' channels. For the complete sample, there is no correlation between the gas-line strengths or their 
ratios and either the sources' dust continuum flux or infrared slope. Modeling shows that fractional abundances of 10"'-10"" can 
explain the line intensities observed, although significant modeling uncertainties remain. For DO Tau, which we model in great detail, 
we find that two very different temperature and density disk structures produce very similar lines for the same underlying abundances. 
Instead, it is the gas properties (particularly and ^gas) and the projected kinematics (determined by M^, sin i) that have the largest 
impact on line appearance; and these system parameters are not well-constrained by current dust models, but instead will be probed 
directly with ALMA. 

Conclusions. Unresolved observations of the dust continuum provide neither a unique nor a complete picture of protoplanetary disks. 
Instead, gas-line measurements and resolved observations of dust and gas alike are needed to arrive at a full picture. 

> ■ 

a^ 

1 . Introduction ular gas content that is left unprotected by the dust is subject to 

^—1 increased photoionization and photodissociation processes. 

• Many young pre-main-sequence (PMS) stars of a few solar One key development in the evolution of disks is the growth 
^ [ masses or less (T Tauri stars) are surrounded by di sks of gas of dust particles from sub-micron sizes first to larger sizes, and 
and dust in the m idst of forming planetary systems ("B eckwithl ultimately planets. The SEDs of many disks exhibit clear evi- 
ll996tlDutrev et al .(,2007). Observations with new instruments in denc e of this growth to m i cron and even ce ntimeter sizes (e.g. 
the sub-millimeter and infrared wavelength regimes have pro- Bec kwith & Sargenill99lt iTesti et aT]l2003h . When destructive 
vided a wealth of information about the dust content in these collisions are negligible, this grain growth process inherently de- 
disks. Physical models are now fitted successfully to the ob- pletes the small-grain population, which is the dominant source 
served spectral energy distributions (SEDs), reveali ng details on of UV opacity in the disk. In addition, as the dust particles grow, 
their large-scale radial and vertical structure (e.g. ICalvet et al.l they are e xpected to decouple from the gas and settle to the disk 
l2002h . Meanwhile, the millimeter spectral slope and the silicate midplane (IWeidenschillingll997l:lDullemond & Dominikl2004l) . 
emission features at 10 and lOfim indicate the amount of small- leaving the gas more exposed to both X-ray and UV irradiation 
scale dust processing, including grain growth (e.g. Natta et al. in the upper disk layers. In a recent anal ysis of S pitze r IRS data 
l2007t and references therein). These dust properties determine for a large sample of T Tauri stars. iFurlan et al.l ([20061) find ev- 
the overall opacity of the disk. A broad range in properties, in idence of this sedimentation of dust particles to the midplane. 
turn, implies different self-shielding abilities from the strong X- They conclude from the infrared spectra that the small dust con- 
ray and ultra-violet (UV) radiation fields emitted by the cen- tent has already decreased by 2-3 orders of magnitude in the 
tral star. As this high-energy radiation penetrates the disk more surface layers of many disks in the Taurus star-forming region, 
or less efficiently, possibly irradiating deeper disk layers or a Since the processes of grain growth and dust removal (via set- 
greater proportion of the available gas reservoir, then any molec- tling) both correspond to a lower disk opacity, we can now test 
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the hnk between the dust properties and the photoprocessing ef- 
fects on the molecular gas composition. 

Observations of molecular lines are required to probe 
the chemical content and evolution of the molecular gas. 
To date, only about three dozen disks h ave been detected 



1995; Mannings & Sargent 


1997', '2000'; "Dutreyetal 


' '1997; 


Guilloteau et al. 




1999; Simon etal. 2000; Thi et al. 


2001; 



van Kemt?en et al. 2007; Pietu et al. 2007; Chapillon et al. 2008 
Schrever et al.li2008i: iSchaefer et al.n2009i) . and only a handful 
of disks hav e been searched for rarer species like HCO^, HCN, 
and CN (e.g Dutrev et al. 1997; Kessler-Silacci 2004; T hi etaP 
[2004; Kastner et al. 1997, 2008a b). The latter three molecules 
are indicative of the rates at which X-ray ionization and UV pho- 
todissociation is occurring in the intermediate layers of the disk. 

In those observational studies, HCO^ and CN often repre- 
sent two of the strongest emission lines detected toward disks 
around PMS stars. Models predict that the abundances of HCO^, 
HCN, and CN should increase similarly as the X-ray ioniza- 
tion r ate is increased dLepp & Dalgarnolll996l: iGlassgold et al.l 
l2004l) . However, the UV field then further enhances CN, a pho- 
todissociation product, a t the expense of HC N, its predecessor 
dLepp & Dalgarnolll996t lAikawa etaT]|2002h . In this way, we 
can use measurements of HCO^ as a normalization factor for the 
total ionization (or overall irradiation) effect on the disk, while 
the ratio of CN to HCN probes the UV photodissociation pro- 
cess, which is most sensitive to grain growth and dust settling. 

Young forming stars a re the main co ntributors to the UV 
interstellar radiation field (iHerbig & Goo drich 1986). In star- 
forming regions permeated by moderate UV radiation fields 
from newly forming stars, the ratio of CN to HCN is observed 
to be very strongl y enhanced with respect to typical molec- 
ular clou d values ("Fuente et al."1993^, 'HogerheiideetalJ [19951 
iGreaves & Church 1996; Terzieva & Herbst 1998). The same 
shift is predicted by chemical disk models where the ratio of 
CN over HCN cha nges radically in fav or of CN as the stellar 
UV flux increases dvan Zadelhoff" et al.l 2003). Young stars also 
possess strong X-ray activity levels dGudell2008 ). Therefore, the 
brightest HCO^ lines and the largest CN to HCN ratios are ex- 
pected from lu minous sources and optically thin disks. 

Previously, iKastner et alJ d2b08b) summarized the HCO^, 
HCN, and CN line observations available in the literature for 
disks around PMS stars, showing how the chemical impact of 
the X-ray and UV fields is significant in the general disk popu- 
lation and that the line ratios do exhibit tentative trends that cor- 
relate with the expected photoprocessing effects. However, also 
evident is the small number of widely-differing systems that have 
been studied in any detail at all. This makes it difficult to distin- 
guish evolutionary differences from source-to-source variations 
in mass, luminosity, environment, and orientation. Naturally, 
larger homogeneous surveys are neccessary to identify gas-rich 
disks and to provide observational tests of the processes driving 
their physical and chemical evolution. 

Here we present a study of HCO^, HCN, and CN line emis- 
sion toward all the classical T Tauri stars in Taurus with a com- 
pact 1.3 mm continuum flux >75mJy, or equivalently a disk 
mass > 0.014 M0, using the standard flux relation and adopting 
the same mean disk parameters as iDutrev et al.l dJ996). These 
parameters include a disk temperature of 15 K, a dust absorption 
coefficient Ki ^^^ of 0.02cm^g"', and a distance to Taurus of 
140 pc. Consequently, our selection criterion favors warmer and 
more massive disks in the star-forming region. 

The line data were obtained from new single-dish observa- 
tions, as well as taken from the literature, where available. In 



Sect. [2] we describe the details of the observations and in Sect. [3] 
we report the line intensities and detection limits obtained from 
our measurements and from the literature. Sect. [4| investigates 
observational trends between the dust emission properties of the 
disk sample, such as millimeter flux and infrared slope, and the 
gas-line strengths or upper limits. Because none of these tracers 
spatially resolve the disks, and because line strengths depend 
both on the total amount of molecules and on the excitation 
conditions. Sect. [5] employs a modeling approach to constrain 
the HCO^, HCN, and CN abundances. This is done using two 
generic models (Sect. l5T2l i and, for one source, DG Tau, in partic- 
ular, a more detailed approach was taken (Sect.[53J. Sect.[6]dis- 
cusses our results and the limitations of our methods and Sect. [7] 
summarizes the main findings. 



2. Observations 

2. 1 . The Taurus disk sample 

The sample consists of the classical T Tauri stars in the Taurus 
star-forming region that were included in the Spitzer IRS clas- 
sification work by iFurlan et al.l d2006l) and possess a 1 .3 mm 
continuum flux larger than 75 mJv i Beckwith et al.l[T990i) . The 
sources cover a broad range in spectral indices between 6, 13, 
and 25;um, and the shape and strength of their 10-/im silicate 
emission features. The difference in the mid-infrared spectral 
indices (defined as An = ni3-25 - "6-13) serves as an indica- 
tor of dust settling in the inner disk (<1 AU, with smaller An 
suggesting more settling), and the strength of the lO/zm silicate 
feature decreases with a depleting population of small grains 
CFurlan et al.. .2006 ). This characterization of the current dust 
properties in the upper surface layers of the disk allows us to ex- 
plicitly contrast the presence of small grains - the main source of 
UV opacity - with the line strengths and ratios of the molecular 
species that trace UV photodissociation. 

For better odds of a high signal-to-noise gas-line detection, 
we imposed the 75 mJy continuum flux minimum cutoff. This 
criterion results from the millimeter flux being related to the dust 
mass, and the idea that a larger dust mass is indicative of a po- 
tentially large gas reservoir Our source list includes some of the 
brightest young stellar objects (YSOs) in Taurus at millimeter 
wavelengths, many of which already possess HCO^, HCN, and 
CN rotational line detections. Consequently, the cutoff value was 
chosen to be representative of the lowest circumstellar dust mass 
found in the literature for a disk that had been detected in at least 
one of the three molecules. Therefore, our study is not necessar- 
ily probing a lower dust-mass regime, but rather contributes to 
completing the sample of sources down to this flux cutoff, and 
expanding the subgroup of disks with similar dust masses that 
can be compared. 

We restricted the sample to one star-forming region, so that 
our results are relatively insensitive to the initial cloud condi- 
tions. In addition, due to their shared environment, we can as- 
sume that the local interstellar radiation field is similar for all 
disks, and more importantly, that it plays a less significant role 
in irradiating the disk in the UV regime. It is generally assumed 
that the external contribution to UV photodissociation of the up- 
per and outer disk layers is roughly 3-4 magnitudes lower in 
intensity than the stellar excess UV spectru m for T Tauri stars, 
based on observations of five protostars by iHerbig & Goodrich 
(119861) . of which four were located in Taurus. 

The full sample consists of 21 protoplanetary disks, in- 
cluding 13 that had not been observed previously in any of 
our molecular tracers. Complementery molecular line observa- 
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Table 1. Taurus disk sample. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Source 


RA 


DEC 


F] .3 mm 


An 


lOyum 




Sp. 


Hq-EW 






Name 




(J2000) 


[mJy] 




Strength [Mq] 


Type 


[A] 




nO^^" erp s^'l 






Our suhsample with new 


molecular line observations presented here 






CVJ Tan 
VV 1 all 


OzL- 1 A- 1 7 on 


+Zo. lU. J / .o 


OA -u C 


n OA 
U.Vo 


U.lu 


1 HA 




1 


1 .OJ 




CYTau 


04- 17-^^ 77 


+28:20:46.8 


133±11 


0.50 


0.12 


0.55 


M1.5 


/u 


47 




FT Tan 


UH-.Z, J . jy. 1 O 


+24:56:14.3 


130 ±14 


0.31 


0.43 




M3 




U.Jo 


f) 7Q0 


DGTau 


(\A-11-C\A 70 


+26:06:16.2 


443 ± 20 


0.44 


0.00 


0.67 


K6 


i 13 






IQTau 


04:29:51.56 


+26:06:44.8 


87± 11 


0.40 


0.35 


0.68 


M0.5 


o 
6 


0.65 


0.416 


V806 Tau 


04:32:15.41 


+24:28:60.0 


124± 13 


1.46 


0.35 




MO 


OS 


1.30 




UZTau 


04-32-42 96 


+25:52.31.1 


172± 15 


0.58 


0.45 


0.18 


Ml-2 


oo 
8Z 


1.06 


0.890 


CI Tau 


04:33:52.01 


+22:50:30.1 


190± 17 


0.93 


0.54 


0.77 


K7 




0.87 


0.195 


DNTau 


04-35-27 37 


+24:14:59.6 


84 ± 13 


0.33 


0.22 


0.70 


MO 


Iz 


0.91 


1.155 


DO Tau 


04-38-28 58 


+26:10:49.9 


136±11 


0.64 


0.18 


0.70 


MO 


100 


1.20 


< 0.270 


GO Tau 


04-43 -03 08 


+25:20:18.6 


83 ±12 


0.97 


0.40 


0.71 


MO 


81 


0.27 


0.249 


DQTau 


04-46-53 06 


+ 17:00:00.2 


91±9 


0.19 


0.14 


0.43 


MO 


113 


0.91 




DRTau 


04:47:06.22 


+ 16:58:42.9 


159± 11 


0.61 


0.27 


0.38 


K5 


87 


3.00 








Sources with molecular line data already in 


the literature {see 


Table\4i 






V892Tau 


04:18:40.62 


+28:19:15.5 


289 ±13 


-1.22 


0.65 




B9 


13 


20.90 


8.576 


RYTau 


04:21:57.41 


+28:26:35.6 


229 ± 17 


0.46 


1.36 


1.63 


KO 


20 


7.60 


5.520 


TTau 


04:21:59.43 


+ 19:32:06.4 


280 ±9 


1.19 


-0.10 




K7 


60 


15.50 


8.048 


GGTau 


04:32:30.41 


+ 17:31:40.2 


593 ±53 


0.74 


0.64 


1.55 


K7 


54 


1.55 


0.160 


DLTau 


04:33:39.06 


+25:20:38.2 


230 ± 14 


0.24 


0.09 


0.70 


K7 


105 


0.70 


< 1.000 


DM Tau 


04:33:48.73 


+ 18:10:10.0 


109 ±13 


1.93 


0.96 


0.25 


M1.5 


139 


0.25 


< 0.457 


A A Tau 


04:34:55.42 


+24:28:53.2 


88±9 


0.50 


0.35 


0.78 


K7 


37 


0.74 


0.690 


GMAur 


04:38:21.34 


+26:09:13.7 


253 ± 12 


3.19 


1.19 


1.00 


K7 


96 


1.00 


1.241 



Notes. (Cols. 1-3) Source name and 2MASS coordinates; (Col. 4) 1.3 mm continuum fluxes from iBeckwithI 119961) : (Cols. 5-6) Dust settling parameters An and 
the lO^im silicate strengfli from jpurlan et al. ( 2006); (Col. 7) Stellar masses from White & Ghez ( 2001) or Hartiga n et al] <1995h : (Col. 8) Spectral types from 
iLuhman et all )2010l) or lRebuU et al. t2010); (Col. 9) Ha e quivalent widths from Herb ig & Bell 1 1988); (Col. 10) Bolomet ric luminosities fro m Palla & StahieJ 
t2002h . except for DQ Tau and DR Tau, which come from lKenvon & HartmannI 119951) : (Col. 11) X-ray luminosities from lGiidel et alj j2007l) or lDamiani et alj 

(UU). 



tions for those 13 sources are presented here. Table [T| lists their 
2MASS source coordinates (used for our JCMT observations), 
1.3 mm continuum fluxes (Beckwith et al. 1990), mid-infrared 
dust settling parameters (Furlan et al. 2006!), as well as several 
stellar properties that are closely tied to the radiation field. 

2.2. JCMT single-dish millimeter observations 

Using the RxA receiver and the ACSIS backend on the James 
Clerk Maxwell Telescope (JCMT) on Mauna Kea, we observed 
the HCO+ J = 3-2 (267.558 GHz), HCN J = 3-2 (265.886GHz), 
and CN 7 = 2-1 hyperfine line (centered at 226.800 GHz) emis- 
sion toward the subsample of 13 sources. The observations 
were made between September 2007 and December 2OO90lThe 
HCO^ observations were carried out in position-switched mode 
with a reference position at a 900" offset in azimuth. Although 
a 900" position switch is unlikely to be off" the Taurus molecu- 
lar clouds, the observed lines trace relatively dense gas and little 
emission is expected away from the source. For the less abundant 
HCN and CN species, we used the more time-eflicient beam- 
switching mode with a reference position at a throw of 180". 

All observations had a bandwidth of 250 MHz, giving a 
channel resolution of 30.5 kHz (or ~0.034 km s ' for HCO^ and 
HCN, and ~0.040kms-' for CN). The atmospheric opacity was 
between 0.05 and 0.12 during the observations, with typical on- 
source integration times of 20 minutes (aside from one extra- 
deep 60-min observation of HCN toward the brightest source, 
DG Tau, in order to establish a deeper upper limit for the sub- 
sample). 



' JCMT observing programs M07BN06 and M08BN01 (PI: Salter). 



The JCMT beam is 19" at 265 GHz and 22" at 226 GHz, 
much larger than any protoplanetary disk at a distance of 140 pc 
to the Taurus star-forming region ( K enyon et al., 1994). Since 
many young stars are still situated near their parental cloud, 
any line detection needs to be checked for possible c ontamina- 
tion fr om t he extended cloud. Follow ing the strategy o f iThi et al.1 
('2001) and v an Kempen et a l. (2007), every strong detection of 
HCO^ was followed up with observations 30" east and 30" north 
of the object. 

2.3. Data reduction 

The data were calibrated by the RxA pipeline and reduced with 
the Starlink KAPPA software. The edges of each spectrum were 
removed and a low-order (p = 2) polynomial baseline was sub- 
tracted. Multiple integrations on individual lines were averaged 
at this stage. The resulting spectrum was smoothed using a mov- 
ing square 'averaging' box of width 3 channels. We then binned 
the smoothed data by 5 channels. The final data resolution is 
152.5 kHz (or0.17kms-i forHCO+ and HCN, and 0.20 km s ' 
for CN). Finally, the reduced spectra were converted from 
to main beam temperatures T^b using a beam efficiency factor 
of 77 = 0.650 The final spectra are shown in Figs. [T]-[3] and the 
resulting rms noise levels are reported in Tables |2] and |3] The 
analysis of the reduced data was performed using the GILDAS 
CLASS software to take advantage of the available hyperfine 
line-fitting algorithms for blended lines. Results of the Gaussian 
fits to the data, including the peak temperature (r,nb), line widths 
(FWHM), and the integrated line intensities ( jT^ah Sv), are also 

- http://docs.jach.hawaii.edu/JCMT/HET/GUIDE/het_guide/ 
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Fig. 1. HCO^ (J - 3-2) observations of our 13-disk subsample in Taurus. A vertical line above the spectrum indicates the central 
radial velocity ulsr for the detected emission, as determined from a Gaussian fit to the emission profile. In the case of the non- 
detection toward DR Tau, we indicate the expected ulsr based on our fit to its CN emission profile in Fig. [3] 



' ' ' I ' ' ' I ' ' ' I 




■U|_SR (km s '') 

Fig. 2. HCN (7 = 3-2) observations of our subsample in Taurus. A vertical line above the spectrum indicates the central velocity 
found from the detected HCO^ emission in Fig.[Tl except for DRTau where the detected CN emission velocity is used instead. We 
note that HCN emission is not detected toward any of the sources, including the most massive disk DG Tau. 
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0.3 



CW Tau 



0.0 



0.3 



0.0 



CY Tau 




DN Tau 



0.3 
0.0 




0.3 
0.0 

0.3 



I ' ' ' I 
DO Tau 



DQ Tau 



DR Tau 



FT Tau 



GO Tau 



.0 ^^^^"'■^^^^"^^"^^S^ 



0.3 
0.0 




10 Tau 



UZ Tau 



V806 Tau 



T^LSR (km s '') 

Fig. 3. CN(7 = 2-1) hyperfine line observations of each of our 13 disks in Taurus. The three vertical lines above each spectrum 
indicate the expected positions of the three strongest hyperfine lines, located at 226874. 1660, 226874.7450, and 226875.8970GHz. 
The spectra have their velocity axis defined such that the strongest hyperfine component, at 226.8747450GHz, should appear at the 
source velocity (vlsr)- In cases where the source velocity is unknown, the Taurus cloud velocity of 6.0km s"' is indicated instead, 
as a reference. We report CN detections toward CI Tau, DG Tau, DR Tau, and GO Tau only. 



reported for HCO+ and HCN in Table|2]and for the CN hyperfine 
lines in Table[3] 



3. Results 

In Fig. [U we see that only 7 of our 13 sources show de- 
tected emission Unes of HCO+ {J = 3-2) : CI Tau, CW Tau, 
DG Tau, DN Tau, DO Tau, GO Tau, and V806 Tau. The de- 
tected lines have intensities of r^b = \22-A9\ mK and are well 
fit by single Gaussian profiles with widths of ~1 km s ' FWHM. 
Remarkably, no strong double-peaked line profiles characteris- 
tic of inclined disks are found. Only DG Tau and GO Tau show 
wider lines of 1.44 and 1.93kms"\ respectively. These sources 
also represent two of the strongest lines, with respective peaks 
Tmb of 49 1 and 25 1 mK. GO Tau is the only source for which the 
line profile is not well described by a single Gaussian. Instead, it 
consists of a superposition of three, narrow (~0.5kms"') com- 
ponents. Greaves (2004) and Thi et al. (2004) also comment on 
the complex line profiles of GO Tau, attributing some of the 
components to surrounding cloud material. Sect. 15.41 discusses 
the literature and these detections further. 

Fig. H shows the HCO^ spectra at 30" east and 30" north 
of the four sources with the brightest detected HCO^ lines: CW 
Tau, DG Tau, GO Tau, and V806 Tau. In all four cases, hints 
of extended cloud emission are seen. However, only toward CW 



Tau are they of comparable strength and shape as the on-source 
emission, and we no longer consider this a detection of emission 
from the disk. This brings our total disk detections in HCO^ 
down to 6. For the other sources, the ofF-source emission is 
< 30% of the on-source line strength, which we determine is 
sufficiently small to be inconsequential to our trends analysis 
(Sect.lU; while our more detailed modeling of DG Tau considers 
the contributions that remnant envelopes and other circumstellar 
material can add to the line emission (Sect.|5]l. Finally, because 
the remaining HCO^ detections (CI Tau, DN Tau, and DO Tau) 
are without accompanying off'-source measurements and they 
exhibit relatively narrow lines, their disk origin must also be con- 
firmed. In the meantime, because their line widths are broader 
than what is typically expected from extended cloud material 
(~0.3 km s"' for a cloud with a T of 25 K and an H2 number den- 
sity of 10"^ cm"^), we report these sources as detected, in follow- 
ing suit with previous single-dish studies. However, our results 
highlight the relevance and necessity of performing off-source 
measurements during single-dish surveys, even for molecular 
species that preferentially trace much denser material. 

No HCN (3-2) emission was detected toward any of our tar- 
gets, down to 3o- upper limits of 41-1 10 mK in 0.17kms"' bins 
(see Table |2]i. Following the recipe of J0rgensen et al. (2004)), 
this corresponds to upper limits on the integrated intensities 
ranging from 61 to 163mKkm s Clearly, toward the sources 
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Fig. 4. From left to right, the HCO+ {J = 3-2) on-source and 30" offset observations for CW Tau, DG Tau, GO Tau, and V806 Tau. 
These 4 sources show the strongest HCO^ emission profiles in our subsample. A dotted line marks the central radial velocity i/lsr 
in km s ' for a Gaussian fit to the HCO^ on-source emission. In all cases, HCO^ emission is detected at the off positions. However, 
with the exception of CW Tau, it is always at a separate velocity or a weaker intensity, which is why we attribute the emission 
toward CW Tau to the surrounding cloud and not the disk. 

Table 2. HCO^ and HCN line measurements. 



Source 


[km s-'] 


0"rms 

[mK] 


max 
[mKf 


S/N 


HCO+ (3-2) 
peak r,^b FWHM 
[mK]* [kms-i] 


S Tnib Av 

[mKkms-']' 


/ T'mb Sv 
[mKkms-i]'' 


[mK] 


HCN (3-2) 
/ T'mb Sv 
[mKkms-']'' 


CI Tau 


6.56 


57 


181 


3.2 


161 


0.68 


105 


117 ±32 


76 


< 112 


CWTau 


6.83^ 


55 


241 


4.4 


259 


0.39 


92 


108 ± 22 


76 


< 112 


off 30" N 


6.95 


61 


279 


4.5 


270 


0.38 


95 


109 ± 27 






off 30" E 


6.93 


61 


248 


4.1 


251 


0.38 


88 


101 ±24 






CYTau 




62 








1.00 




<92 


110 


< 163 


DGTau 


6.47 


100 


491 


4.9 


424 


1.44 


695 


652 ± 72 


41 


<61 


off 30" N 


6.07 


58 


203 


3.5 


174 


0.43 


89 


79 ±23 






off 30" E 


6.13 


65 


251 


3.9 


247 


0.47 


107 


124 ±29 






DNTau 


8.00 


48 


122 


2.6 


128 


0.44 


56 


60 ± 18 


85 


< 126 


DO Tau 


6.16 


48 


182 


3.8 


127 


0.86 


115 


116±26 


67 


<99 


DQTau 




59 








1.00 




<88 


81 


< 120 


DRTau 




65 








1.00 




<97 


93 


< 138 


FT Tau 




54 








1.00 




<80 


87 


< 129 


GO Tau 


5.05 


47 


251 


5.3 


186 « 


1.93 


327 


383 ±39 


71 


<105 


off 30" N 


7.37 


68 


148 


2.2 


279 


0.17 


46 


51 ±19 






off 30" E 


7.30 


71 


126 


1.8 


155 


1.03 


172 


170 ±44 






IQTau 




80 








1.00 




< 119 


91 


<135 


UZTau 




66 








1.00 




<98 


80 


< 119 


V806 Tau 


5.40 


77 


345 


4.2 


271 


0.97 


297 


279 ± 45 


100 


< 148 


off 30" N 


5.80 


61 


162 


2.6 


201 


0.25 


82 


52 ±20 






off 30" E 


5.59 


68 


164 


2.4 


149 


0.39 


76 


62 ±28 






GO Tau 1'' 
GO Tau 2 
GO Tau 3 


4.40 
5.40 
6.34 


47 
47 
47 


268 
251 
143 


5.7 
5.3 
3.0 


237 
232 
125 


0.54 
0.75 
0.35 


100 
179 
50 


136 ±28 
185 ±28 
46 ± 19 





Footnotes, (a) The maximum temperature value in the binned data; (b) The maximum determined by a Gaussian fit to the fine using the CLASS Gauss fitting 
routine; (c) A sum of the individual channels from the biimed data to be compared to the integrated fit; (d) For non-detections we report a 3<t upper hmit where 
(T= 1.2iTjnis VAViSv and we assume a fine width AV of l.Okms"'; (e) Upper Hmits are reported in the same way as for HCO^; (f) We have determined that the 
emission detected here is not associated with the circumstellar disk; (g) Statistics are reported for a single Gaussian fit, which may not be reahstic for this source; (h) 
Statistics for the individual peaks toward GO Tau, measured as they appear from left to right in Fig.[T] 
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Table 3. CN line measurements. 



Source 


CN (2 3/2,3/2- 


-1 5/2,5/2) 


CN (2 3/2,5/2- 


-1 5/2,7/2) 


CN (2 3/2,1/2- 


-1 5/2,3/2) 


Average Stats 




!^LSR 


Peak 


/ 7'mb Sv 


!^LSR 


Peak 


/ 7'rab Sv 


!^LSR 


Peak 


/ Tn,b Sv 


FWHM 


0"rms 




[kms-'] 


[mK] 


[mKkms"'] 


[kms-'] 


[mK] 


[mKkms"'] 


[kms-'] 


[mK] 


[mKkm s"'] 


[kms-'] 


[mK] 


ClTau 


7.33 


230 


52 ±19 


6.51 


76 


16±19 


5.28 


65 


58 ±35 


0.42 


63 


CWTau 






<94 






<94 






<94 


1.00 


58 


CYTau 






<97 






<97 






<97 


1.00 


60 


DGTau 


7.63 


198 


171 ±94 


6.45 


279 


414 ± 124 


5.03 


294 


76 ±30 


0.82 


63 


DNTau 






< 104 






< 104 






< 104 


1.00 


64 


DOTau 






<90 






<90 






<90 


1.00 


56 


DQTau 






<98 






<98 






<98 


1.00 


61 


DRTau 


9.65 


119 


56 ±26 


11.06 


297 


87 ±20 


11.93 


110 


24 ± 16 


0.31 


55 


FTTau 






<88 






<88 






<88 


1.00 


55 


GOTau 


2.56 


55 


101±41 


4.07 


103 


37 ± 19 


5.24 


116 


129 ±32 


1.03 


42 


IQTau 






<91 






<91 






<91 


1.00 


57 


UZTau 






<79 






<79 






<79 


1.00 


49 


V806 Tau 






< 102 






< 102 






< 102 


1.00 


63 



Notes. The velocities reported here are for a mlsr = centered at the theoretical frequency for the strongest hyperfine component, at 226.8747450 GHz. All reported 
peaks were determined from a Gaussian fit. All upper limits are given by 3 X 1.2 (Tnns VAV^v following jorgensen et at.. (.2004.) where AV is set to 1.0km s"' . 



Table 4. IVIolecular line fluxes. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Source 


FlIco 


'^CO 


Ref. 


Fhco+ 


HCO+ 


Ref. 


Fhcn 


HCN 


Ref. 


FcN 


CN 


Ref. 


Name 


[Jy kms"'] 


Trans. 




[Jy kms"'] 


Trans. 




[Jy kms"'] 


Trans. 




[Jy kms"'] 


Trans. 





AA Tau 


12.4 


2-1 


2 


1.9 


1-0 


8 


1.8 


1 


-0 


8 


3.8 


1-0 


8 


CI Tau 








2.4 ±0.7 


3-2 


1 


<1.6 


3 


-2 


1 


2.6± 1.5 


2-1 


1 


CWTau 


8.3 


2-1 


2 


<1.8 


3-2 


1 


<1.6 


3 


-2 


1 


<1.9 


2-1 


1 


CYTau 








< 1.9 


3-2 


1 


<2.3 


3 


-2 


1 


<2.0 


2-1 


1 


DGTau 


40.0 


1-0 


3 


13.5 ± 1.5 


3-2 


1 


<0.8 


3 


-2 


1 


13.7±5.1 


2-1 


1 


DLTau 


<2.1 


3-2 


2 






















DM Tau 


6.5 


2-1 


4 


4.1 ±0.5 


3-2 


5 


2.0 ±0.3 


1 


-0 


6 


8.7 ±0.4 


2-1 


6 


DNTau 


14.5 


2-1 


2 


1.2 ±0.4 


3-2 


1 


<1.8 


3 


-2 


1 


<2.2 


2-1 


1 


DOTau 


31.0 


2-1 


2 


2.4 ±0.5 


3-2 


1 


< 1.4 


3 


-2 


1 


<1.9 


2-1 


1 


DQTau 








<1.8 


3-2 


1 


< 1.7 


3 


-2 


1 


<2.0 


2-1 


1 


DRTau 


4.3 ± 1.2 


3-2 


5 


3.7 ±0.4 


4-3 


9 


<1.9 


3 


-2 


1 


3.5± 1.3 


2-1 


1 


FTTau 








< 1.7 


3-2 


1 


<1.8 


3 


-2 


1 


<1.8 


2-1 


1 


GGTau 


5.8 ±0.2 


2-1 


6 


6.7 ±0.9 


3-2 


6 


2.8 ±0.7 


3 


-2 


6 


8.4 ±0.4 


2-1 


6 


GM Aur 


9.3 ±2.7 


3-2 


5 


9.9 ±1.2 


4-3 


9 
















GOTau 


4.1 ±1.2 


3-2 


5 


7.9 ±0.8 


3-2 


1 


<1.5 


3 


-2 


1 


5.5±1.9 


2-1 


1 


IQTau 








<2.5 


3-2 


1 


<1.9 


3 


-2 


1 


<1.9 


2-1 


1 


RYTau 


4.3 ± 1.2 


3-2 


5 


1.7 ±0.7 


4-3 


9 


1.6 


1 


-0 


8 


0.8 


1-0 


8 


TTau 


3.0 ±0.2 


2-1 


7 


3.1±0.1 


1-0 


7 


1.9±0.8 


2 


-1 


10 








UZTau 








<2.0 


3-2 


1 


< 1.7 


3 


-2 


1 


<1.6 


2-1 


1 


V806 Tau 








5.8 ±0.9 


3-2 


1 


<2.1 


3 


-2 


1 


<2.1 


2-1 


1 


V892 Tau 
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in our Taurus subsample, flie HCN lines are at least a factor of 3 
weaker than the HCO^ lines (and a factor of 12 weaker toward 
DG Tau where a much deeper integration resulted in a lower 
noise level). 

Two sources, DG Tau and GO Tau, show clear detec- 
tions of the CN(2-1) triplet at frequencies of 226.8741660, 
226.8747450, and 226.8758970 GHzQ Two additional sources, 
CI Tau and DR Tau, show tentative detections. No emission was 
detected from the other targets where the typical rms noise level 
was 57 mK in 0.20km s"' bins. Two suggestive peaks in the DO 
Tau spectrum at the correct position for the CN hyperfine com- 
ponents are likely noise. 



' The JPL Molecular Spectroscopy Catalog http://spec.jpl.nasa. gov| 



4. Trends 



The previous section shows 6 detections of HCO^ emission, 
no detections of HCN emission, and 4 detections of CN emis- 
sion from our subsample of 13 observed sources. Looking at the 
complete sample of 21 sources yields the following statistics for 
detections of HCO+ (14/21), HCN (5/21), and CN(8/21) toward 
the T Tauri disk population thus far in Taurus. Half of the re- 
ported detections are taken from the literature (see Table |4|i and 
these often include rotational transitions other than the ones ob- 
served by us. However, in the present section we focus on gen- 
eral trends in the detected emission strengths, and we assume 
that the different transitions correspond to line strength differ- 
ences of no more than a factor of ~2. 
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Fig. 5. A log plot of the integrated line intensities in Jy km s"' versus several disk parameters that probe the total dust mass (Fumm), 
dust settling via the mid-infrared spectral slope (An), amount of small grains via the 10/zm silicate emission feature strength (Si- 
strength), and the total gas mass (inferred from the '^CO line flux). Blue squares represent our JCMT observations, red circles 
represent data from the literature, and black arrows are used to indicate all upper limits. A green star symbol represents T Tau, 
whose line emission is contaminated by a large remnant envelope in the single-dish data plotted. The upper-right panel illustrates 
the dynamic range of the available gas reservoirs. (Color plots are available online) 



Our goal is to investigate whether the line strengths (and the 
upper limits) from the Taurus sample show any correlation with 
disk or stellar properties. For this purpose, we converted the line 
strengths to integrated line fluxes (in Jykms"^) in order to re- 
move the dependences on beam size. We used the formula: 



• 7.354 X 10" ■ rMB[Jy/K] 



(1) 



where D. is the primary beam in arc seconds, A is the observed 
wavelength in cm, and Tmb is the observed main beam tempera- 
ture in K. The assembled line fluxes for the entire sample are 
provided in Table |4] To trace the disk properties we use the 
millimeter continuum flux /^i.3m m, the difference in the mid- 
infrared slope An as defined bv iFurlan et alJ (20061). and th e 
10-;um silicate feature strength, also from iFurlan et al ] (I2006h . 
Respectively, these trace the disk (dust) mass, dust settling, and 
small-dust content. Stellar properties are probed by stellar mass 
M*, Ha equivalent width, bolometric luminosity Lboi, and X-ray 
luminosity Lx- We stress that all quantities are observed quan- 
tities, and that a relation (or lack thereof) between the HCO^ 
line strength and the 1.3 mm continuum flux, for example, does 
not immediately imply a connection (or lack thereof) between 
HCO^ abundance and disk (dust) mass. To draw the latter con- 



clusion, detailed modeling of the emission mechanisms is re- 
quired, which is the topic of Sect. |5] 

One of the guiding questions in our research was, whether 
disks with a lower dust content or with more dust settling have 
larger gas volumes subject to UV irradiation, and therefore show 
comparatively stronger lines of HCO^ (more ionization) and CN 
(more photodissociation). In the following figures, we will ex- 
plicitely investigate these relations. Fig. |5] plots the observed 
HCO+, HCN, CN, and '^CO line fluxes against Fum^, An, 
lO-jum siUcate feature strength, and '^CO line flux. The latter 
quantity serves as a tracer of the total molecular gas content. We 
choose to plot the logarithm of these quantities because we are 
primarily interested in large-scale trends and because we want 
to suppress differences of a factor of ~2 due to variations in 
source distance, inclination, luminosity, and plotted line transi- 
tions. No immediate trends are present in the panels of Fig. |5] 
Instead, we find that: (1) the millimeter continuum flux alone 
is not a clear predictor of line strength since both lower- and 
higher-mass disks exhibit detections and non-detections alike 
in all species; (2) disks with and without significant dust set- 
tling exhibit equally strong line strengths all around; and (3) a 
larger '^CO reservoir is not necessarily accompanied by stronger 
HCO+, HCN, and CN fines. These results hold for a '^CO line 
flux spanning a factor >10 in dynamic range. 
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Fig. 6. A plot of the integrated line intensities in Jy km s ' versus several stellar parameters from the literature, including: the stellar 
mass (M*), the He equivalent width (an accretion tracer), the bolometric luminosity (Lboi), and the X-ray luminosity (Lx)- See 
Fig.|5]for an explanation of the symbols. (Color plots are available online) 



We then considered the same HCO^, HCN, and CN line 
strengths versus j^n,. An, lO-yum silicate feature strength, and 
'^CO line flux, but this time only after normalizing the line fluxes 
to the '^CO line flux (not shown). This relation represents the 
line flux per unit gas mass, if we assume that the '^CO line fluxes 
are reliable tracers of the total amount of molecular gas in each 
disk. No obvious trends were visible there either. Finally, we ex- 
plored the equivalent plots of line flux, but then normalizing to 
^^i.3mm, or the line flux per unit dust emission. Again, no trends 
emerged. Together, the results show that there is no clear corre- 
lation between the observed line fluxes of HCO^, HCN, or CN, 
and the disk dust mass, the degree of settling, the amount of 
small particles in the disk, or the disk gas mass. Echoing our 
earlier remarks, we emphasize that only detailed modeling of 
the emission lines can prove the presence or absence of such a 
correlation; but judging from the figures that plot general trends 
only, no strong correlation is expected from this sample. 

If the molecular line flux does not depend strongly on the 
properties of the disk, perhaps it depends on the properties of 
the star In Fig.|6]we plot the line fluxes and upper limits against 
stellar mass. Ha equivalent width, bolometric luminosity, and 
X-ray luminosity. Similarly, we also explored the results of nor- 
malizing the line fluxes with respect to the '^CO line flux and 
^i.3mm- While the UV flux could be expected to depend strongly 
on the mass accretion rate, no homogeneous set of mass accre- 
tion rates is available in the literature for this sample, and these 



can be expected to be variable as well. Instead, we use the Ha 
equivalent width as a tracer and find no correlation. Indeed, no 
trends are apparent in any of the plots, suggesting that also such 
a simple connection between the stellar radiation and the emer- 
gent HCO+, HCN, or CN line flux is absent. 

In Fig. |7] we plot the specific relations that served as the 
primary motivation for this study: the ratio of CN over HCN 
line fluxes and the HCO^/'^CO flux, respectively tracing the de- 
gree of photodissociation and the degree of ionization, versus the 
gas-to-dust ratio (represented by '^CO/Fumm) and the chang- 
ing mid-infrared slope An, representing dust settling. As was the 
case for our other plots, no significant correlations are apparent. 
We note that the decreasing line strength with a larger gas-to- 
dust ratio in the upper-left panel of Fig. |7] is the likely effect of 
the incorporation of the '^CO factor into each ratio. 

In summary, we find no correlation between the HCO^, 
HCN, and CN line fluxes (or their ratios) and any tracer of the 
disk properties or those of the stellar radiation field. The line 
flux does not seem to be affected by any of the investigated 
parameters. We conclude that the d etails of the input r adiation 
field, such as UV and Lya strengths (iBergin et al.ll2003]) . may be 
the deciding factor in the resulting line fluxes. Other contribut- 
ing factors include the inner versus outer disk contributions and 
the temperature and density structure in the line-emitting region, 
both of which are discussed in Sect. |5] 
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Fig. 7. Upper panels: The de- 
gree of ionization (traced by 
HCO+/'^CO) versus the gas-to-dust 
ratio ('"'CO/Fi 3mm), the difference 
in the mid-infrared spectral slope 
between 6, 10, and 25//m (An), and 
the strength of the 10//m silicate 
emission feature. Lower panels: 
Similar plots for the photodissoci- 
ation effect (traced by CN/HCN). 
See Fig. |5]for an explanation of the 
symbols. (Color plots are available 
online) 



5. Modeling the molecular emission 

5.1. Disk models 

The previous section investigates the observational correlation 
between the measured HCO^, HCN, and CN line intensities and 
upper limits, and the disk dust observables, such as millimeter 
continuum flux and infrared slope. While the emergent line in- 
tensity depends on the underlying molecular abundance, other 
factors including the disk density and temperature structure af- 
fect the emerging lines through molecular excitation and line ra- 
diative transfer This section addresses how the modeled abun- 
dances of HCO^, HCN, and CN (that can explain the observed 
emission) are related to the disk dust properties. Rather than de- 
veloping an ab initio description of the disk structure and asso- 
ciated molecular chemistry, this section employs two indepen- 
dent models obta ined from the literature ( Robitaille et al. 200^ 
llsella et ani2009l) as starting points, and calculates the molecu- 
lar abundances (assumed constant throughout the disk except in 
regions of freeze-out) that reproduce our line observations. 

In the first modeling approach, we make use of the online 
SED fitting tooQ of .Robitaille et al.. (.2006.) , the best-fit param- 



eters of Robitaill e et al. I (I2007h and the visual extinction val- 



from White &Ghez ('2001)0 The continuum radiative trans- 
fer code of Whitnev et al. (2003) produced the two-dimensional 
density and temperature structure for the best-fitting model for 
each source. In some cases (i.e. CI Tau, DO Tau, DR Tau, and 
FT Tau), the Robitaille models include remnant envelopes; these 
are not plotted in the figures but are included in the line calcula- 
tions. However, since their H2 number densities are < 10^ cm"^, 
the envelopes are not expected to contrib ute significantly to the 
HCO ^, HCN, or CN line emission (e.g. Hogerheiide & Sandelil 
l2000h . This first method relies entirely on spatially unresolved 
continuum data, keeping all other stellar and disk parameters 
free, even in cases where these properties m ay be well known. 

Our second approach uses a model by llsella et aP (l2009h . 
which explicitely takes into account the spatial distribution of 
the millimeter continuum emission observed with the CARMA 



Online SED fitting tool, 'http://caravan. astro.wisc.edu/protostarsl 
' No visual extinction was included in iRobitaille et alj i2007h and 
therefore our best-fit model may differ slightly from the one listed there. 



interferometer. These authors approximate the vertical temper- 
ature structure of the disk with the two-layer description by 
IChiang & Goldreich ( 1997), and fit the grain size and opacity to 
resolved 1 .3 mm data. In our adaptation of their models, we omit 
the hot surface layer because it contains insignificant amounts 
of molecular gas, and we extrac t the disk's interior temperature 
from Fig. 7 o f Isella et al. I (l2009l) . The surface density is obtained 
from their Eq. 9 and Table 5. We then calculate the local hydro- 
dynamic equi librium scale height following Eqs. 3 and 4 from 
iHughes et al.l (|2008|). We truncate th e models at the transition 
radius, defined by Isella et aTl (l2009l) as the location where the 
disk becomes optically thin to the stellar radiation (Sect. l53] dis- 
cusses the effect of extending the disk further). 

Figs. [8] and |9] show the resulting temperature and density 
structures for our sources. For some sources, only one type 
of model is shown b ecau se of the availab ility of models in 
IRobitaille et"al] (|2007|) and llsella etal] (l2009b . Table |5] summa- 
rizes the parameter fits for both models. As can be immediately 
seen for the four sources for which both a Robitaille and an 
Isella model are available (i.e. CY Tau, DG Tau, DN Tau, and 
DR Tau), widely different disk descriptions apparently provide 
equally good fits to the continuum data. 

The vertical height of the Robitaille models is often 4-6 
times smaller than the Isella models. The radial extents of the 
disks are comparable for both approaches to within a factor of 
<2. The resulting masses differ by factors of 2-10, with the 
Isella models always producing more massive disks. The tem- 
perature profiles of the Robtaille models are more detailed (by 
definition) while the Isella models at large radii are close to 
isothermal. In general, the temperature between the two mod- 
els differs by a factor of 3-5 in the midplane. In addition to these 
parameters, other factors that will affect the emergent lines are 
the disk inclinations and stellar masses, both of which influence 
the widths of the lines. 

To calculate the resulting molecular line emission, we pop- 
ulate each model disk with gas using the standard 100:1 gas-to- 
dust ratio, a mean molecular weight of 4.008 x 10"^^ g, and con- 
stant H2 relative abundances of 1 x 10 for HCO+, 2 x 10 " for 
HCN, and 1 x 10"^ for CN, based on the fractiona l values in the 
warm mol ecular layer for the theoretic al models of lAikawa et al.l 
(I2OO2I) and lvan Zadelhoff et all (l2003h . In regions where the dust 
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Table 5. Summary of the stellar and disk properties for the best-fitting dust models for our 13 disks. 



Source 






Robitaille Model 








Isella Model 








Av 




Md 






Fl 3mm 




M,( 




i 


1 3mm 




[mag] 


[Mq] 


no"^ Mr,^ 




[°] 


[mjy] 


[Mq] 


no"^ Moi 




[°] 


[mJy] 


ClTau 


1.80 (1.80) 


0.4 


1.9 


91 


57 


105 












CW Tau 
























CYTau 


0.05 (0.03) 


0.9 


1.4 


97 


32 


23 


0.4 


6.92 


197 


51 


117 ±20 


DGTau 


1.60 (1.41) 


1.5 


3.9 


158 


50 


254 


0.3 


41.7 


89 


18 


317 ±28 


DNTau 


0.49 (0.49) 


0.6 


1.0 


92 


32 


39 


0.4 


1.86 


125 


30 


93 ±8 


DO Tau 


4.88 (2.23) 


0.8 


3.2 


104 


57 


103 












DQTau 
























DRTau 


0.51 (0.51) 


0.8 


3.2 


104 


18 


103 


0.4 


6.31 


86 


37 


109 ± 11 


FT Tau 


0.00 ( - ) 


0.2 


0.8 


62 


50 


36 












GO Tau 














0.6 


7.10 


670 


25 


57 ±8 


IQTau 


1.44(1.44) 


1.0 


3.0 


101 


63 


74 












UZTau 














0.3 


4.79 


260 


43 


126 ± 12 


V806 Tau 

























Notes. For the RobitaiUe models, the visual interstellar extinction Ay magnitudes are for the best model fit, only using the measured values from IWhite & Ghe3 
i2001h as an input parameter to the SED fitting tool (in parenthesis). The My, values listed for both models ai'e the derived values for each method. For the Isella 
models, the outer dust radius is from Table 5 in'I sella et all 12009'). defined as where the disk becomes optically thin to the stellar radiation. For UZ Tau, the Isella 
model in fact concerns only the spectroscopic binary component UZ Tau E. The 1.3 mm continuum fluxes (f i.3mm) listed for our best-fitting Robitaille model are 
the SED values of the fit; and for the Isella models we list the resolved (interferometric) dust continuum measurement. 



temperature (assumed equal to the gas temperature) drops be- 
low the CO ice evaporation temperature of 20 K, these abun- 
dances are reduced by a factor 10^ to account for CO freeze-out. 
When the H2 number density drops below 10^ cm"^, the molecu- 
lar abundances are set to zero, effectively defining the disk edge. 

The gas kinematics follow a cylindrical Keplerian velocity 
field with stellar masses corresponding to each model fit, as in- 
dicated in Table |5] A Doppler broadening factor of O.lSkms ' 
is also factored into the line calculations. The statistical equilib- 
rium molecular excitation a nd line radiative transfer was solve d 
using the RATRAN code dHogerheiide & van der l2000h . 
and the emerging line emission was convolved with the appropri- 
ate Gaussian beams. The resulting line profiles are plotted along- 
side the disk models in Figs. [8] and |9] 

5.2. Comparison of fixed-abundance models 

As is immediately obvious from Figs.|8]and|9] the models have 
difficulty reproducing the detected emission lines. In four cases 
(CI Tau, DG Tau, DN Tau, and DO Tau) the models predict lines 
that are weaker than observed; and in two cases (GO Tau and UZ 
Tau), the models predict lines that are stronger than the detected 
line (for GO Tau) or the obtained upper limits (for UZ Tau). In 
the remaining four cases (CY Tau, DR Tau, FT Tau, and IQ Tau), 
the predicted lines are consistent with the achieved upper limits, 
although for CY Tau the Isella model produces an HCO^ line 
that violates the upper limit. 

Where both a Robitaille model and an Isella model are avail- 
able, the Isella model always produces lines that are stronger 
and narrower than those of the Robitaille models. The larger 
stellar masses, by factors 2-5, of the latter, and the larger disk 
masses of the former, by factors 2-10, contribute to this differ- 
ence. Inclination also plays a significant role, with more face- 
on orientations leading to stronger lines (cf. DO Tau and DR 
Tau, which are fit with the same Robitaille disk structure, but 
have respective inclinations of / = 57° and / = 18°; DR Tau has 
predicted lines stronger by a factor of ~3). Interestingly, of the 
six sources from our sample modeled by Isella et al. (2009), the 
three with detected line emission have i < 30° (DG Tau, DN 
Tau, and GO Tau), while the sources with / > 30° (CY Tau, DR 



Tau, and UZ Tau) are undetected (all have similar M*). Perhaps 
the narrower line profiles have helped to make these sources de- 
tectable. 

Given the general mismatch between the predicted line in- 
tensities and widths, and the observations, it is not possible to 
draw conclusions about the HCO^, HCN, or CN abundances 
this way. Simply scaling up or down the abundance will not 
result in a match (to the line shape); only for GO Tau do the 
abundances appear to lie within a factor of a few above the true 
values. Furthermore, as illustrated by the case of DR Tau, two 
different disk models (but with more comparable Mj, Rd, and 
^^i.3mm values) produce very similar lines. These degeneracies 
make it difficult to derive reliable conclusions about the HCO^, 
HCN, and CN abundances. Instead, more detailed modeling of 
individual sources may be required. 

5.3. The specific case of DG Tau 

Of all the sources in our sample, DG Tau offers the best case 
to obtain a detailed model. Its 1.3 mm continuum is brighter 
by a factor >2.5 than any of the other sources and it emits 
the strongest HCO^ and CN Unes. An extensive literature 
also exists on DG Tau describing sub-m illimeter single-dish 
(Schuster et al. 1993; Mitchell et al."1994|) a nd interferometer 
data (Kitamura et al. 1996ab; Dutrev et aP 119961: iTesti et all 
Isella et al T l2009r Its milUmeter continuum emission is 
compact with ~80% originating from within 95 AU dlsella et al.l 
2009). The best-fit Robitaille model obtained in the previous sec- 
tion (see Table |5] Robitaille database model number 3017659) 
matches other literature estimates of the DG Tau disk parameters 
(Table©. Only the disk mass (Md = 0.042 Mq) for the Robitaille 
model is significantly larger than the literature values (0.015- 
0.025 Mq). In contrast, the Isella model offers a disk mass value 
lOx greater (0.4 Mq). The Isella model also overestimates the 
stellar mass (M* = 1.5 Mq) with respect to most literature values 
where M* ranges from 0.3 to 0.8 M©; the Robitaille model yields 
0.3 Mq in comparison. To model DG Tau we settle on a central 
stellar mass of 0.8 Mq, which is on the high end of the literature 
values for this object but provides the best fit to the line profiles, 
as discussed below. 
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Figure [10] reproduces the '^ CO(6-5), '^CO(l-O) and 
C"^0(2-l) lin e observations of 'Schus ter et aP jl993h and 
iKitamura et alj ([l996a). The bright CO isotopologues sug- 
gest a significant gas reservoir. Interferometric imaging of the 
'^CO(l-O) and (2-1) lines in the lit erature reveal a gaseo us 
disk structure o f 600 AU in extent (IKitamura et all Il996allbl: 
iTesti et aT]l2002h . about 4-6 times the size inferred from dust 
emission. Accordingly, we adopt 600 AU for the outer radius of 
the gas disk, extrapolating the initial Robitaille mode l outwards. 
In add ition, the '^CO(2-l) emission observed bv iTesti et al.l 
(I2OO2I) is consistent with Keplerian rotation around a star of 
0.67 ± 0.25 Mq, oriented perpendicular (within ±15°) to the 
highly collimated jet sys tem, which is inclined 38° with re- 
spect to the line of sight dEisloffel & Mun dt 1998). Therefore, 
we adopt an inclination of 25°, which also gives the best fit to 
the C"^0 (2-1) line profile. We already show in Fig. [10] — and we 
will discuss later in this section — the minimal effect from inter- 
vening cloud or remnant envelope material in this line. 

With the stellar mass (M* = 0.8 Mq), inclination (/ = 25°), 
and outer gas radius (R^ = 600 AU) now fixed, we use a con- 
stant fractional CO abundance of 2 x 10 and C'^O abundance 
of 4x10"^ (except when r<20K, where an abundance 10^ 
times lower is used) to calculate the simulated line profiles with 
RATRAN using the extended and modified Robitaille model. To 
fit the lines, particularly C'^O, we find that we need to increase 
the gas temperature by a factor 1 .7, suggesting that the line emis- 
sion originates in layers where Tg^s > Tdust- We adopt this same 
scaling for the gas temperature for all species, but neglect its 
effect on the scale height. 

For these model parameters, the HCO^ line can be repro- 
duced very well for a disk-averaged abundance of 2 x 10 " with 
respect to H2. This is lower than theoretica l predictions for the 
war m emission layer in many T Tauri disks dAikawa et al.ll2002L 
Ivan Zadelhoff et al. 20031), but not unlike the beam-averaged 
fractional abundances (~10"''-10"'^) inferred from observa- 
tions of_disks around several high-mass (Herbig Ae/Be) stars 
rsee lThietal]l200l . A low mean abundance is also especially 
surprising given DG Tau's powerful jets, which emit significant 
X-ray radiation (Giidel et al. 2008). Apparently, the disk is suffi- 
ciently shielded to retain a low ionization degree. An upper limit 
for the HCN abundance of 5 x 10 is found, while a value of 
8 X 10"^° is obtained for CN. This high CN/HCN ratio of > 160 
suggests efficient HCN dissociation in the bulk of the disk, which 
is more consistent with its mid-infrared characterization indicat- 
ing fewer small grains and some dust settling. Here we note the 
importance of using abundance ratios rather than line intensity 
ratios (which contain opacity and excitation effects): the ratio of 
the integrated intensities of CN/HCN used in Sect. [3] is > 12, a 
factor of 1 3 smaller than the underlying abundance ratio found 
here, albeit consistent. 

The disk abundances above have been derived using the 
density and temperature structure of the best-fitting Robitaille 
model, with noted modifications. The results are summarized in 
Table [7] Substituting the description of Isella instead (but keep- 
ing the same values for M* and /), the emergent line intensities 
are lower by up to a factor of ~16 for HCO^ for the same abun- 
dances, and by smaller factors of ~5 and ~10 for HCN and CN, 
respectively. The large line intensity differences result also from 
the differences in Rg.^s for each model. Whereas the Robitaille 
power-law model for the temperature and density description 
(and sharp outer dust edge) leads easily to extrapolation to larger 
radii, the Isella model, with its exponentially tapered density 
structure, does not. Thus, to compare the two disk structures 
in a uniform way, we plot the temperature and density profiles 



for each model for the inner 160 AU only in Fig. [TT] Then, in 
Fig- Owe plot the predicted lines. Unlike Figs. [8] and [9] we now 
assume identical gas kinematic properties (M*, /, and Rg^s)- The 
predicted lines for each model in Fig.[T2]are now strikingly simi- 
lar The differences in the two emerging line intensities are lower 
by up to a factor of 2 for HCO^ and much closer for HCN and 
CN. This suggests that the large CN/HCN ratio of > 160 found 
for the underlying abundances is independent of the temperature 
and density details of the adopted model (for the inner 160 AU); 
and leaves only the missing details of the input radiation field 
uninvestigated. 

Finally, DG Tau is not an isolated sourc e. The environment 
aroun d the star is dominated by an optical jet dEisloffel & Mund^ 
Il998l) . a strong molecular outflow (Mitchell et al. 1994), an ex- 
panding circumstellar envelope ( Kitamura et al.ll 1996 b). and in- 
tervening cloud material (this work). The result of this con- 
fused environment is most clearly evident when compar i ng the 
'^C O (3-2) line observati ons presented in Schuster et alJ d 19931) 
and [Mitchell et alJ d 19941) . which exhibit equally bright line in- 
tensities and significant wings at the on-source position and three 
separate offset positions. We chose to omit the '^CO observa- 
tions from the fits in Fig. [TO] since they distract from the disk 
emission. However, we do determine that the '^CO lines are 
about 3x stronger than predicted for our disk model, suggest- 
ing contributions from a surrounding cloud with a CO column 
density of A^co ~ 6 x 10'^ cm"^ and line width of 0.3 km s"' for 
an adopted cloud temperature of 25 K and H2 number density of 
10** cm"-' (or A^co ~ 1 x lO'* cm"^ for an H2 number density of 
10^ cm"^). For these typical cloud densities, we do not expect 
significant HCO+ (3-2), HCN (3-2), or CN(2-1) emission (or 
absorption); and both Fi g s . [4] and [TOl confirm this. 

To determine those line contributions from the cloud, plotted 
in Fig. [To] as well as the cloud fractional abundances and column 



Table 6. DG Tau star and disk properties from the literature to 
compare with the two disk models. 



Property 


Isella 


Robitaille 


Literature 


References 




Model 


Model 


Range 




Sp. Type 


MO 




K1-K7 


2, 10, 12, 23 


M, [Mo] 


0.3 


1.48 


0.3-0.8 


4, 5, 17, 20 


Age [Myr] 


0.1 


1.89 


0.3-2.0 


5, 11, 13 


Teff [K] 


3890 


4560 


3890-4395 


5, 7, 11, 17 


[Lo] 


1.70 




1.07-3.62 


1,4, 5, 18 


R* [Rg] 


2.87 


2.427 


2.13-2.8 


3, 11, 17 


M [Moyr-'] 


4.1e-^ 


4.5e-^ 


1.2-20e-' 


11, 17, 18 


Md [10-2 Mo] 


41.7 


3.9 


1.5-2.51 


2, 6, 9, 13 


Rdu., [AU] 


89 


158 


80-300 


9, 16, 19 


^gas [AU] 


160 


600 


600 


16, 24 


ni 


18 


25 


18-90 


3, 9, 14, 16 


Av [mag] 




1.6 


1.41-1.6 


10, 22 


Dist. [pc] 


140 


140 


128-156 


10, 21 



References. (1) 'Akeson et alJ (2OO50; (2) I Andrews & Williams? (200^; 
(3) Appenzeller et al. (2005); (4) Basri et al. 1 1991); (5) Beckwith et al. ( 199(|); 
(6) Beckwith & Sargent 119911); (7) iBouvier et aU >1995i) ; (8) Cieza et al] 
(2005); (9) Dutrev et al. ( 1 991); (10) iFurlan et alj f2006l); (11) iHartigan et al] 
119951); (12) IHessman & GuentheJ fT997h; (13) Honda etal. (200^; 
(14)llsellaetn] <2009l) ; fl5) lKenvon & Hartmannl {l995); (16) Kitamura et al] 
)1996bl); (17) 'Mohanty et al. (2()0 |); ( 18) Muzerolle et al. (200^; 
(19) lRodm aniiet al. (2006); ( 20) Ta m uraet alj <1999.); (21),yink p yic.&JurkB 
( 2007); (22) White & Gheg ilMIh ; (23) IWhite & Hillenbrand jI5^ ; and 
f24) iTesti et aH2003.) . 
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Table 7. Determined abundances for our best-fit DG Tau model. 



Molecule DG Tau Theoretical 



Disk Fractional Abundances (w.r.t. Hi)" 



"CO 


2.0 X lO^" 


l.Ox lO-"* 


'^co 


3.3 X 10-* 


1.7 X 10-" 


c"o 


4.0 X 10-^ 


2.0 X 10-' 


HCO+ 


2.0x 10^" 


1-100 X 10-" 


HCN 


5.0x 10^" 


1-100 X 10-" 


CN 


8.0x 10-'" 


1-100 X 10-" 


Cloud Fractional Abundances (w.rt. H2)'' 


"CO 


2.0 X 10-* 


8.0 X 10-^ 


"CO 


3.3 X 10-'' 




c"o 


4.0 X 10-^ 




HCO+ 


8.0 X 10-' 


8.0 X 10-' 


HCN 


8.0 X 10-' 


4-20 X 10-' 


CN'^ 




3-30 X 10-' 


Cloud Column Densities (cm ^) 


"CO 


6.0x 10'« 




"CO 


l.Ox 10'^ 




c"o 


1.2 X lO"* 




HCO+ 


2.4 X 10" 




HCN 


2.4x 10" 




CN 







Notes, (a) The DG Tau abundances are constant throughout the disk (or disk- 
averaged), whereas the theoretical disk abundances from Aikawa et al. (2002) 
and , van Zadelhoff et all ^20031) represent the ranges expected in the warm molec- 
ular layers only; (b) The theoretical cloud values are from iTerzieva & Herbsll 
(1993) : (c) The onhne RADEX program does not yet include CN (to calculate 
the cloud contributions). 

densities listed in Table |7] we used the RADEX0 online calcu- 
lator. We modeled the cloud as a cold intervening layer mov- 
ing with a radial velocity of 6.1kms-^. We note that in many 
of the literature observations, strong absorption is seen near a 
velocity of 5.8-6.2 km s-', with several studies reporting these 
values as the source velocity. We emphasize here that our obser- 
vations of HCO^, with a critical density 3 orders of magnitudes 
larger than the CO observations in the literature, are a much bet- 
ter tracer of the disk content, establishing the source radial ve- 
locity at 6.47 km s-'. In addition, we confirm that the emerging 
HCO^ line predicted for our intervening cloud fits the detections 
observed at both DG Tau offset positions (~0.25 K peak centered 
at 6.1 km s-', refer back to Fig.|4]i. The effect of the intervening 
material on the observed '^CO, HCO^, and HCN lines is shown 
in Fig.[TOl while '^^CO and C'**0 exhibit no differences, and the 
CN cloud predictions are not yet available in the RADEX pro- 
gram. 

5.4. Notes on individual sources 

V806 Tau, also called Haro 6-13, has the second strongest 
HCO+ line after DG Tau, but is undetected in HCN and CN. 
It is a single MO sta r (Rebull et a l. 201 0; Luhmanet al. 2010). 
With an Ay = 1 1 .2 (jpurlan et al.l l2006h . its optical extinction 
is much larger than our other sources, which might help ex- 
plain its u nique silicat e em ission feature and positive spec- 
tral slope. iFurlan et al.l (12006) comment that the silicate fea- 

* RADEX is an online, one-dimen s iona l, non-LTE radiative transfer 
code developed by van der Tak et alJ ( 12007 ) to calculate the molecular 
line intensity, opacity, and excitation temperature. For more informa- 
tion, see http://www.sron.rug.nl/~vdtak/radex/radex.phpi 



tures are reminsicent of transitional di sks, although it also 
possesses a high mass accretion rate (IWhite & Hillenbrand 
12004). Some extended HCO^ emission is apparent in Fig. |4] 
However, our single-dish radial velocity of 5.40 km s-' is con- 
sistent with the value of 5. 10 km s-' from CO interferometric 
observations (Schaefer et al. 200 ^. V806 Tau's low disk m ass 
of 0.01 Mo (Andrews & Wilhams 2005: Honda et al.l l2006h and 
large 400 AU radius (Robitaille et al. 2006) suggest that much of 
the disk of V806 Tau should be UV illuminated, in contrast to 
what our observed upper limit for CN would suggest. 

GO Tau shows complex line profil es, with peaks at a wlsr of 
4.4, 5.1,5.5, and 6.3 km s-'. lThi et alj (1200 1) report '^CO emis- 
sion peaks at 5.2, 5.5, 6.2, and 7.1 km s'\ and '^CO emission 
at 4.3 and 7.0 km s-'. They attribute the 5.5 and 6.2 km s-' 
components to surrounding cloud emission. However, we as- 
sign the emission between 4-6 km s-' to the disk of GO Tau , as 
suggested by the interferometric observations of iSchaefer et al.l 
( 12009.) a nd .Andrews C2007,). Fo r GO Tau, the accretion rate is 
very low (iHartmann et al.ll 19981) . and the amount of dust settling 
is small if we draw upon its mid-infrared spectral slope, suggest- 
ing lower rates of UV photodissociation and ionization should be 
occurring. However, it exhibits one of the brightest HCO^ lines, 
in stark contrast to the rest of the sample since its fainter 1.3 mm 
continuum flux just straddles our cutoff and yet the source also 
appears to be associated with a large, and dense, gas reservoir 

DR Tau also has a complicated circumstellar environment, 
with '^CO emission lines at 6.8, 9.1, 10.0, and 10. 3 km s^' 
and ' -^CO emission at 6.9 km s-' and near 1 1 kms-' (|Thi et al.l 
l2001h . SMA interferometric observations by lAndrewsl (|2007|) 
also show strong emission centered on 10.5 km s-', suggesting 
that this is the correct source velocity and lending credence to 
our assertion that the CN emission likely originates in the disk. 
Our CN spectrum shows a triple-peaked 5.4cr feature centered 
at 11.2 km s-' (Fig. [3]i. Our model results suggest that 'stan- 
dard' HCO^ and HCN abundances are consistent with the non- 
detections of the lines. The CN detection in the absence of the 
other lines, on the other hand, indicates a significant CN en- 
hancement. 

CW Tau is clearly surrounded by dense cloud material, as 
witnessed by the equally strong HCO^ emission on- and off- 
source. Our HCO^ observations (see Fig. HJi illustrate how, in a 
crowded star-forming region, measurements at offset positions 
can be both relevant and useful even for molecular species that 
preferentially trace much denser material. 

CY Tau, DQ Tau, IQ Tau, and UZ Tau do not show 
gas emission lines in our data. Interestingly, they span the full 
dust classification and morphological sequence of iFurlan et al.l 
(2006) with CY Tau and DQ Tau showing rather flat, decreas- 
ing mid-infrared SEDs, and IQ Tau and UZ Tau showing evi- 
dence for a small grain population. The mass accretion rates - 
which contribute to the stellar UV excess - range from very low 
( I O-** Mo yr- ' ) for C Y Tau to av erage ( 1 0-^ Mq yr- ' ) for DQ Tau 
and UZ Tau ( Giidel et al.l |2007'). However, both DQ Tau and UZ 
Tau E are spectroscopic binaries that exhibit pulsed accretion 
events on periods of weeks (Bas ri et al.lll997h . and their higher 
reported accretion rates may overestimate the average, quiescent 
values. 

6. Discussion 

To return to previous work, iKastner et al.l (|2008b") showed plots 
of HCO^, HCN, and CN line ratios for several PMS stars, show- 
ing a tentative correlation between the HCO^, HCN, CN, and 
'^CO line ratios. In Fig. [13] we reproduce their plots and their 
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data points (without error bars), and add our own line ratios. 
Overall, we find that the trends persist: CN is typically stronger 
than HCN, CN is also stronger than HCO^, and the photodis- 
sociation rate (as probed by CN/HCN) is roughly constant re- 
gardless of the HCN relative line strength. Our contributed line 
ratios consist largely of upper and lower limits and therefore do 
not specifically challenge or confirm the trends by probing dif- 
ferent regions in the plots. While this complete PMS sample in 
Fig. [T3] includes sources covering a range in age, mass, and ra- 
diation field, the line strengths of HCO+, HCN, and CN rela- 
tive to one another do reveal the importance of the ongoing UV 
photodissociation and X-ray ionization processes in these disks. 
The results, however, are still limited by small numbers statis- 
tics, numerous upper limits, and different rotational transitions 
(that may trace different regions of the disk). 

The motivation of this study was to determine whether disks 
with a higher degree of dust settling, or with a decreased dust 
content, have higher abundances of CN and HCO^ reflecting 
larger degrees of photodissociation and photoionization. Our 
data are inconclusive. In Sect. |4] we found that the HCO^, HCN, 
and CN line fluxes (or their ratios) do not depend on any other 
disk or stellar parameter such as millimeter flux, infrared slope, 
silicate feature strength, stellar spectral type, etc. Sect. |5] shows 
that detailed SED-based models have intrinsic degeneracies that 
preclude straight-forward modeling. And even a detailed model, 
tailored to the case of DG Tau, does not provide unambiguous 
estimates of the HCO^, HCN, and CN abundances. 

This suggests two possible ways forward. In the first, spa- 
tially resolved observations of both the dust continuum and the 
line emission can be used to obtain in situ measurements of the 
molecular abundances. ALMA will be a powerful instrument for 
an analysis like this. By addressing localized disk regions, rather 
than the emission integrated over the entire face of the disk, more 
simpUfied modeling approaches can be used (not unlike what 
is presently state-of-the-art analyses for photon dominated re- 
gions, PDRs). In addition, spatially resolved observations pro- 
vide many more constraints on the underlying disk structure, 
such as the extent and surface density. This approach addresses 
the question of how molecular line emission and underlying disk 
structure are interrelated. 

The second approach does not focus on the details of the 
underlying disk structure, but rather the details of the star's ra- 
diation field, which irradiates the disk. Accurate determination 
of the spectral type and luminosity, the UV and X-ray emis- 
sion characteristics, and their time dependences should lead to a 
better understanding of the response of the molecular gas reser- 
voir to the incident radiation. High-resolution, high signal-to- 
noise observations are time-consuming but necessary for bet- 
ter stellar characterization of PMS stars, while veiling and 
other circumstellar en vironmental effects provide challenges that 
have been overcome (iMerin et al ] |2004H h erczeg & Hillenbrand' 
I2OO&) . Recent studies have even argued that the continuum UV 
spectrum for T Tauri stars must be much weaker, due to the total 
fraction of the stellar FUV flux that is emitted in the Lya line 
alone. Since this fraction can range from 30% up to 85% (for 
TW Hya), Lya provides an additional source of photodissoci- 
ating power that va ries from source-to -source and is extremely 
difficult to measure dBergin et al.ll2003l) . Which of these two ap- 
proaches is most fruitful will, of course, depend on which fac- 
tor dominates the molecular line emission: the underlying disk 
structure or the stellar irradiation. Since these may be interre- 
lated, both approaches may prove necessary. 
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Fig. 13. Integrated intensity line ratios for HCO^, HCN, CN, and 
'^CO plotted alongside other disk systems around PMS stars in 
the literature. Circles represent data points extracted from Fig. 3 
of (Kastner et al. 2008b), squares represent data points from our 
complete Taurus sample, and black arrows represent upper/lower 
limits. 



7. Summary 

We surveyed 13 classical T Tauri stars in low-7 transitions of 
HCO^, HCN, and CN to compare the gas structure and chemi- 
cal abundances within planet-forming disks that possess similar 
dust masses. We found, conversely, a wide variety of molecular 
gas properties. For this sample in the Taurus star-forming region, 
we report 6 new detections of HCO^ (3-2), new detections of 
HCN (3-2), and 4 new detections of CN(7 = 2-1). These data 
double the pool of previously known detections, bringing the to- 
tal detection statistics toward the 21 brightest (at 1.3 mm) disks 
in Taurus to: 14 for HCO+, 5 for HCN, and 8 for CN. 

Overall the HCO+, HCN, and CN line ratios for our Taurus 
disk sample are consistent with the trends identified toward other 
disks around PMS stars found in the literature, as initially plot- 
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ted bv lKastner et alJ (l2008allbl) . In general, CN is more prevalent 
than HCN, which suggests that the bulk of the detected emission 
originates in a UV photodissociating region. Additionally, the 
fractional molecular ionization ratio, as traced in only slightly 
denser regions by the HCO^ line, is also enhanced. Both trends 
agree with the narrow emission lines observed in our sample, 
which trace the outermost regions where the disk is optically 
thin to both the stellar and interstellar radiation fields. 

Despite this disk-to-disk agreement in the line ratios toward 
the general population of unresolved disks, the gas-line proper- 
ties reveal no observed chemical photoprocessing effects due to 
the dust properties or several stellar parameters, which was the 
main motivation for this research. We do not see CN and HCO^ 
enhancements (via brighter lines) in sources whose mid-infrared 
spectral slope and silicate emission features indicate more grain 
growth and dust settling (leading to a lack of UV shielding). In 
addition, stellar parameters like X-ray luminosity do not seem to 
influence the observed line intensities of ionization tracers such 
as HCO+. 

The next step was to derive the underlying molecular abun- 
dances using two dust models in the literature that are then pop- 
ulated with theoretical values for the fractional molecular abun- 
dances. Models for the sample as a whole illustrate the impor- 
tance of the M* sin / and /?gas factors in gas-line abundance stud- 
ies; parameters that are poorly constrained by the dust properties, 
but critical to proper line fits. Along these lines, we found dur- 
ing detailed modeling of the source DG Tau that the underlying 
abundances were less dependent on the temperature and density 
details in the adopted dust models. We conclude that better char- 
acterization of the stellar parameters (M*), the radiation field it- 
self (UV and Lya), and spatially-resolved line observations (Rgas 
and are necessary to constrain the molecular gas content and 
evolution. 

Large sample statistics are still a challenge. It remains a fu- 
ture task for ALMA, with its the resolution and sensitivity ad- 
vancements, to resolve these disks (including inner and outer 
disk differences), their gaseous surface density profiles, and the 
chemical signatures of changes in the inner dust structure. Only 
then will we be able to examine how photodissociation, ioniza- 
tion, and freeze-out processes affect the surface density of the 
gas by comparison to resolved dust observations. 
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Fig. 8. Results of the general analysis. From left to right: the Robitaille disk structure; the line predictions for HCO^, HCN, and CN 
(in red for the Robitaille model, green for the Isella model); and the Isella disk structure (when available). The sources from top to 
bottom are: CI Tau, CY Tau, DG Tau, DN Tau, and DO Tau. For the disk structures, filled contours correspond to the temperature 
profile (in K). The temperature levels are identical for all disks, designated at: 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 75, 100, 
150, and 250 K. A color bar is provided in Fig.|9] The contour lines indicate the log H2 number density (in cm ■^), indicated at whole 
number integrals from 3 (disk surface/edge) to 12 (typical midplane density). Models are only pictured for the disks included in the 
Robitaille and Isella samples. (Color plots are available online) 
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Fig. 9. Continued from Fig. [8] The sources from top to bottom are: DR Tau, GO Tau, FT Tau, GO Tau, and IQ Tau. A color bar 
is provided here for the filled temperature contours and applies to all disks. Models are only pictured for the disks included in the 
Robitaille and Isella samples. 
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Fig. 10 . Line predictions for our best-fitting model of DG Tau. The '^CO (6-5) and C '^O (2-1) data a re taken from lSchuster et aP 
11991 where their T* scale is equivalent to our r^b scale. The '^CO(l-O) is from iKitamura et al] ([T996a), and since we were 
unable to establish a conversion from Jy to K, we compare the line shape only. We have indicated with a thick blue line our fits for a 
600 AU disk at an inclination of 25°. Indicated in red is how the ('^CO, HCO^, and HCN) line profiles are affected by the absorption 
or excess emission from a cold, foreground cloud at a radial velocity of 6. 1 km s ' . No differences are visible for the '^CO or C'^O 
lines. The effect of the foreground cloud on the CN line is not available; but the expected hyperfine line-splitting is calculated. 
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Fig. 11. The radial and vertical structures for the inner 160 AU of each disk model, for direct comparison. The panels at left show 
the Robitaille model (without the extended gas reservoir), and the panels at right show the full Isella model (with exponential taper). 
The upper panels indicate the temperature structures, and the lower panels compare the density structures for each model. (Color 
plots ai^e available online) 
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Fig. 12. Similar to Figs. [8]and|9] the line predictions for the two models for the inner 160 AU, assuming identical abundances, Tkin, 
M*, and /, in order to uniformly compare the predicted emission. A solid blue line represents the Robitaille solution, whilst a red 
line is the Isella prediction. The corresponding temperature and density structures are given in Fig. [TT] In addition, as a dashed 
red line, we plot the most extreme solution to the Isella model, based on the provided error bars and the largest possible fractional 
abundances predicted for disks. For both models, the inner 160AU structure emits only a small fraction of the total emission 
necessary to recreate the observed lines, but the two predictions are very similar. 
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